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Evidence for multiple impurity bands in sodium-doped silicon MOSFETs 
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We report measurements of the temperature dependent conductivity in a silicon MOSFET that 
contains sodium impurities in the oxide layer. We explain the variation of conductivity in terms of 
Coulomb interactions that are partially screened by the proximity of the metal gate. The study of 
the conductivity exponential prefactor and the localisation length as a function of gate voltage have 
allowed us to determine the electronic density of states and has provided arguments for the presence 
of two distinct bands and a soft gap at low temperature. 

PACS numbers: 71.23.Cq, 71.55.Gs, 71.55. Jv, 72.15.Rn, 72.20.Ee, 72.80.Ng, 73.20.At, 73.40.Qv 



Since the invention of the silicon MOSFET, under- 
standing the influence of impurities, especially sodium 
contamination, on device performance has been a prior- 
ity and continues to provide a rich system for investi- 
gation by experimental and theoretical physicists alike. 
The electronic properties of sodium doped MOSFETs 
were first studied by Fowler and Hartstein^in the 1970s. 
They reported a single, broad peak in the subthreshold 
drain current against gate voltage and attributed it to the 
formation of an impurity band induced by the presence 
of sodium ions near the Si-Si02 interface. Further stud- 
ies of narrow channel devices (~ 100 nm) demonstrated a 
series of reproducible sharp peaks^, while later exper- 
iments found evidence for resonant tunneling between 
localised states in the channeli^*^ For sufficiently low 
impurity concentrations, the overlap between neighbour- 
ing localized electron wavefunctions and consequently the 
hybridisation of their excited states is predicted to be 
reduced^, splitting the single impurity band observed at 
high concentrations into the ground and excited bands 
as modeled by Ghazali.4 Increasing the resistivity of the 
silicon substrate reduces the scattering from acceptors 
at the Si-Si02 interface, allowing the possibility for such 
a band splitting to be experimentally observed in the 
transport. In this paper, we will present evidence for the 
observation of two separate impurity bands with a soft 
gap, based on analysis of the temperature dependent con- 
ductivity below 20 K. 

The device we used is a MOSFET fabricated on a 
(100) oriented p-silicon wafer and was subsequently pat- 
terned in the circular Corbino geometry to eliminate Hall 
voltages and possible leakage paths. The effective gate 
channel length and interior width were respectively 1 /im 
and 346 /zm. A high resistivity wafer (lO^fi.cm) pro- 
vided a background concentration of less than 10^^ cm"'^ 
of boron corresponding to a mean distance between im- 
purities of 1 ixm. A 35 nm gate oxide was grown at 950 °C 
in a dry, chlorine-free oxygen atmosphere. The phospho- 
rous implanted and aluminium sputtered contacts were 
highly metallic and Ohmic at all temperatures investi- 
gated. Sodium ions were introduced onto the oxide sur- 
face by immersing the device in a 10~^ N solution of high 
purity sodium chloride (99.999%) in de-ionised water. 
The surface of the chip was dried with nitrogen gas and 



an aluminium gate subsequently evaporated. To observe 
or remove the low temperature conductivity structures, 
the mobile ions are drifted through the oxide to the Si- 
Si02 interface, or returned to the Al-Si02 interface by 
applying either a -1-4 V or a —4 V DC gate-substrate bias 
for lOmin at 65°C before the device is cooled down to 
helium temperature at which sodium looses its diffusiv- 
ity in the oxide. All measurements were performed using 
standard low-noise lockin techniques with an amplifier 
of lO^V/A. The AC excitation was maintained at 15 /j,V 
with a frequency of 11 Hz. Suitable RC filters were em- 
ployed to eliminate any DC offset from the amplifier. The 
gate voltage was controlled by a high resolution digital 
to analog converter. All experiments were performed in 
an "^He cryostat and the temperature was measured by a 
calibrated germanium thermometer. 

Fig. 1 shows the conductivity a of our device at 300 mK 
versus gate voltage for the case where the sodium ions 
had been drifted to the Si-Si02 interface. Two groups of 



peaks appear clustered around Vg 



-2V and -0.5 V 



and separated by a region of low conductivity and lim- 
ited by noise. The origin of the peaks themselves will be 
discussed later. Following a — 4V drift, no structure was 
detectable over the full range of gate voltages investigated 
but a difference in threshold voltage of 0.2 V was found 
at 77 K between the characteristics of the device follow- 
ing -I-4V and — 4V drifts. This was attributed to the 
presence of mobile charges close to the Si-Si02 interface 
at a concentration of 3.7x lO^'^ cm~^ corresponding to a 
mean impurity separation of 16 ± 1 nm. In a reference de- 
vice where no sodium was introduced, no sub-threshold 
conductivity peaks and no shift of the threshold voltage 
appeared for any drift conditions investigated. Tunnel- 
ing through the oxide and other leakage currents were 
discounted as the gate leakage current was below 50 f A 
at 4.2 K and was approximately constant over the range 
of gate voltages used. We remark also that no hysteresis 
and thus no charging effects that have been reported in 
similar devices^ were observed here. 

The presence of two distinct ranges of Vg where 
peaks appear suggests the possibility of a split impu- 
rity band (Fig. 1). Such a splitting into a ground and 
an excited state is expected to happen for low-doping 
concentratior&iS if one takes into account the overlaps 




FIG. 1: The source drain conductivity versus gate voltage at 
300 mK following a +4V drift and a -4V drift (dotted line). 



between impurity wavefunctions and uses a multi-band 
formalismii. The conductivity of a Si-MOSFET is not 
directly related to the density of states and the fact that 
we see two regions of high conductivity separated by a 
region of low conductivity is only an indirect indication 
that the device density of states consists of two bands sep- 
arated by a gap. In fact, through the Kubo formalism, 
conductivity tends to be related to local paths through a 
disordered device but density of states is a global prop- 
erty. In order to show that the density of states splits 
into two bands, we have looked at the temperature de- 
pendence of the conductivity at a series of different gate 
voltages. For all gate voltages studied, the conductivity 
decreases non-monotonically as temperature is lowered 
(Fig. 2). In the range 1 K to 20 K, we observe the charac- 
teristics of hopping conduction so that the conductivity 
ct(T) is fitted to the generalised equation : 



a = aoT 



(1) 



where Tq and ctq depend on gate voltage. 

The best values for 7 and s are found by minimizing 
the value of the reduced chi-square deviation 
standard procedures. Figs 2(a) — (b) show the optimum 
values determined for four gate voltages, one point in 
the threshold region {Vg = 0.1 V), the upper band {Vg = 
-0.24V), the gap {Vg = -1.48V) and the lower band 
[Vg = -2.26 V). The resulting fits for a{T) are valid over 
three orders of magnitude in a (Fig. 3). Studies in gate 
voltage show that 7 = 1.98 ± 0.04 and s = 0.39 ± 0.02 
for Vg below 0.25 V. Above this point, the value of s de- 
creases rapidly towards s = 1/3 for T > 4K and the 
range is better described by Mott hopping conduction^^. 
It is worth noticing that the smallest values of s are ob- 
tained for band centre regions — O.SV^ < Vg < —0.5V 
and -2.2V < Vg < -2.1V for which the hopping 
lengths are relatively smaller and the Coulomb interac- 
tions weaker. 
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FIG. 2: a) Temperature dependence of the conductivity for 
Vg = 0.1 (o), -0.24 (•), -1.48 V (★) and -2.26 (<), b) Varia- 
tion of the reduced with s for 7 = 2 for the gate voltages 
listed showing minima at s = 0.412, 0.406, 0.394 and 0.385 
respectively and c) Variation of the reduced with 7 for 
s equal to the optimum value found in 2(b), for the appro- 
priate gate voltages, displaying consistency with a minimum 
at 7 = 2 (dotted line). Lines in b) represent the exponent 
for the Mott hopping regime (left) and Efros and Schklovskii 
regime (right). 



= 1.83 and s = 0.412 
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FIG. 3: Temperature dependence of the conductivity for 
V"g = 0.1 (o), -0.24 (•), -1.48 V (★) and -2.26 (<) for the op- 
timum values of 7 and s defined in Eq. 1. For clarity, curves 
are shifted downward respectively by 0, 1, 2 and 3 from their 
original values. 



We emphasise that the use of temperature dependent 
exponential prefactors allows the fine distinction between 
Mott (s = 1/3), Efros-Schklovskii (s = l/2)i^andthe 
regime under study (s ^ 0.39). The finding that 7 ~ 2 
is consistent with the formulation for a{T) given by Allen 
and Adkinsi^ii^ if it is rederived for the 2D case. This 
gives the conductivity in terms of both the localization 
length ^ and the density of states at the Fermi level 
n{Ep). Accounting for the fact that Coulomb interac- 
tions between electrons in different localized states mod- 
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ifies the density of states close to the Fermi level so that 
n{E) = Nq \ E - Ef I^, we obtain : 



p+1 

cT = (Ta T~'^(tS) e-(^)^ with 7 = 2 (2) 



^0 



P + 3 
p+1 



(3) 



(4) 




where Aq is a constant depending on the electronic prop- 
erties of bulk silicon and ks the Boltzmann constant. 

The conditions given by Alleni^ for the use of the 
equations Eqs. 2-4 are satisfied in our device. Recent 
calculations showed that sodium ions in the oxide may 
trap either one or two valence electrons against the Si- 
Si02 interface and that the wavefunctions of the localized 
states remain hydrogen-like.^^ Also, hopping conduction 
is present in our device but Coulomb interactions are 
such as the only mobile electrons are found in an energy 
band of few fc^T around the Fermi level. The resulting 
Coulomb gap for s = 0.39 (i.e. n{E) ~ \E- Ep \°-^°) 
is much sharper than for the Efros regime where n{E) ~ 
I E — Ep I . This behaviour has been predicted by Blanter 
and Raikh— while studying 2D systems localized by dis- 
order. They have shown that a metallic gate close to the 
interface provides image potentials that modify the den- 
sity of states to the form n(£:) = n(Ep)+No\E-Ep\'^/^ 
at T = K and thus gives an exponential dependence of 
rp-OA ^]-^g conductivity. This behaviour is explained 
by the fact that the oxide thickness plays the role of the 
screening length and that initial and final states become 
electrostatically independent at low temperature when 
the hopping length R becomes greater than twice the 
oxide thickness d. This then produces a crossover from 
the Efros to the Mott regime. The localization length 
in our device was approximately estimated by fitting the 
conductivity with p = 0, using the Mott formula for 
ksTo and taking the 2D value for the density of states. 
For Vg = —0.4 V, this gives ^ ~22nm and a hopping 
length R ~63nm at IK. Our device has a gate ox- 
ide of d = 35 nm and the Coulomb interactions may 
be screened by the electrostatic gate as i? ~ 2d. This 
demonstrates the device may well be in the regime de- 
scribed by Blanter and Raikh. Fits of the conductivity 
using J — 2 and s — 0.4 are still valid for Vg smaller 
than 0.25V and from IK to 18 K typically, so we will 
proceed with these values. In this regime, the density of 
states at Ep is given by : 

n{Ep) = iio{Ep) '- V (5) 



FIG. 4: Variation of the DOS at Ep with gate voltage. The 
dashed and dotted lines represent respectively the background 
density due to the conduction band tail and the upper band 
tail. The inset is the relative variation of the DOS to the 
background in the same units. 

where is the permittivity of silicon. 

These equations have been derived at T = OK. At 
higher temperature but for fc^T V{d,)^ where V is 
the coulomb potential energy, the gap is partly filled. 
The Coulomb interactions then become negligible for en- 
ergies below V{d) and the density of states saturates such 
as no{Ep) = no{Ep + V{d))ii^ This condition decreases 
the value of v from 1 down to 0.05. Combining (3) and 
(4), the parameters Nq and ^ were expressed in terms 
of To and cfq, values easily accessed experimentally. The 
density of states was extracted using equation (5) and (6). 
To find the value of Aq, we used S = 8.91 eV for the effec- 
tive deformation potential of acoustic phonon&LS- as well 
as 3800 m.s^^ for the speed of surface acoustic phonons 
in siliconSli. This gives Ao = 6.69 x 10~^^ e^./i^^.m^. 

The density of states is shown in Fig. 4. The value 
we find is more than one order of magnitude lower than 
the pure 2D case 1.6 x 10^^ eV~^.cm~^ with valley 
degeneracies). The large background is predominantly 
due to the conduction band tail which spreads over the 
full range of Vg studied. The upper band also has a sig- 
nificant tail. The presence of density of states tails is a 
common occurance in disordered systems with low im- 
purity concentration and localized wavefunctionsSiiSSiSi 
but its linear shape for < F is unusual. It has 
been attributed to the formation of regions of constant 
local potential energy at the Si-Si02 interface and con- 
taining a random number of chargesJ^ Two regions of 
higher density are superimposed on the background den- 
sity and correspond to the upper and lower groups of 
peaks. This confirms that the structure observed in the 
conductivity is due to the presence of two separate bands. 
By numerically subtracting the background density and 
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FIG. 5: Variation of the localisation length (•) and the trans- 
mission coefficient at IK (o) with gate voltage. 

integrating over the appropriate gate voltage and by sup- 
posing a linear relation between the gate voltage and the 
surface potential energy, we estimate the upper band con- 
tains approximately 3 times the number of states as the 
lower band. This value is an upper bound as energeti- 
cally deeper states could not be accessed experimentally. 

The variation of the localization length ^ (Fig. 5) fol- 
lows that of the density of states, showing the two bands 
do correspond to the more conductive regions. The value 
of ^ decreases rapidly when approaching the threshold 
voltage. This is expected as the region between Vg = OV 
and Vg ~ 0.4 V is both in the band tails of the conduc- 
tion band and the upper band, the conduction band edge 
being well above 0.4 V. Thus, this region is still a region 
of strong localization. Nevertheless, the value of ^ is ex- 
pected to rise once the Fermi energy crosses the conduc- 
tion band edge. Finally, taking into account fitting errors 
as well as the discrepancy in the value of the deformation 



potential, we estimate the localization length within 5 % 
and the density of states within 6 %. From the same 
derivation which gave (1), (2) and (3), the transmission 
coefficient is obtained : 

The experimental variation of the transmission coeffi- 
cient in gate voltage (Fig. 5) shows that the conductivity 
of the two bands mostly comes from the higher mobility 
of the states of energies within the upper and lower bands 
and less from the increase in the density of states. 

In conclusion, we have observed an unusual hopping 
regime with an exponent 0.4 which results from the 
screening of the Coulomb interactions by the metal gate. 
We have shown that, consequently, both the localization 
length and the density of states can be extracted from 
the temperature dependence of the conductivity. This 
analysis has given strong evidence for the existence of 
two separate bands and a soft gap at low temperature. 
This may result from the splitting of the impurity band 
into a lower and an upper band in presence of Coulomb 
interactions. The formation of the two bands results it- 
self from the presence of a low concentration of sodium 
impurities close to the Si-Si02 interface. The conditions 
for the observation of such a formation may be the cre- 
ation of deep but well separated impurity potentials at 
the interface, resulting in electron localization. Finally, 
the electron screening may be sufficiently weak and the 
disorder not too strong to allow a Mott-Hubbard transi- 
tion to take place. The two bands could then possibly be 
Hubbard bands. 
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